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ABSTRACT 

The progenitor systems of Type la supernovae (SNe la) are still under debate. 
Based on recent hydrodynamics simulations, non-degenerate companions in the single- 
degenerate scenario (SDS) should survive the supernova impact. One way to distinguish 
between the SDS and the double-degenerate scenario is to search for the post-impact 
remnant stars (PIRSs) in SN la remnants. Using a technique that combines multi- 
dimensional hydrodynamics simulations with one-dimensional stellar evolution simula- 
tions, we have examined the post-impact evolution of helium-rich binary companions 
in the SDS. It is found that these helium-rich PIRSs (He PIRSs) dramatically ex- 
pand and evolve to a luminous phase (L ~ 10 4 Lq) about 10 years after a supernova 
explosion. Subsequently, they contract and evolve to become hot blue-subdwarf-like 
(sdO-like) stars by releasing gravitational energy, persisting as sdO-like stars for several 
million years before evolving to the helium red- giant phase. We therefore predict that 
an sdO-like star should be detectable in the central regions of SN la remnants within 
star-forming regions, if the SN la progenitors evolved via the helium-star channel in the 
SDS. These He PIRSs are predicted to be rapidly rotating (v rot > 50 km s _1 ) and to 
have high spatial velocities (linear ^ 500 km s -1 ). Furthermore, if supernova remnants 
have diffused away and are not recognizable at a later stage, He PIRSs could be an 
additional source of single sdO stars and/or hypervelocity stars. 



Subject headings: binaries: close, — methods: numerical, — stars: evolution, — stars: 
subdwarfs, — supernovae: general 
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INTRODUCTION 



The catastrophic explosions of Type la supernovae (SNe la) are of great importance in probing 
the history of the universe. It is generally believed that SNe la are caused by thermonuclear 
explosions of carbon-oxygen (CO) white dwarfs (WDs) in close binary systems, bu t their intrinsic 
variety a nd nature of their progenitor systems is still under debate (Livio 200ol: Hillebrandt & 



Niemever lfflOol : iRuiz-Lapuentelboid : Iwang fc Hanlboi2h . 



C urrent mainstream progenitor scenarios include the single-degenerate scenario fSDS; W helan 



&Ib 



en 



1973 



Nomotol 1982) and the double-degenerate scenario fDDS; 



Iben fc Tutuko 



1984 



Web- 

bink[l984j). In the SDS, the binary companion is a non-degenerate companion which could include 
stars of many different stel lar types, including main-sequence (MS) sta r s, red giants (RGs), helium 

;_Wang 



(He) stars, and M-dwarfs ( Hachisu et al. , 19991; Ivanova fc Taam 



2009 



2004 



Hachisu et al. 



2008 



20ld : IWang &; Hanll2010al ; IWheelerll2012h . However, these non-degenerate companions 



et al. 

are usually hydrogen - or helium-rich. The observational u pper limit of stripped hydrogen after SN 
impact is < 0.01M Q (lLeonardll2007l : IShappee et al.ll2013bl ). Thus, the absence of the hydrogen that 
should appear in SN la spectra poses a problem for the SDS. The DDS instead results from merg- 
ers of two CO WDs with total mass greater than the Chandrasekhar mass, avoiding the hydrogen 
problem. These violent e yents may lead to acc retion-induced collapse to neutron stars instead of 
thermonuclear explosions (jNomoto fc Ibenlll985l ). However, estimates of the delay time distribution 
(DTD) based on observed supern ova rates are cons i stent with a large fraction of events being due 
to double degenerate progenitors (jMaoz et al.ll201ll . 120121 ) . 



Based on recent hydr odynamics simulations of SN la impact on binary companions in the 
SDS, including grid - based (tMarjetta et ali boool: IPan et al.ll201ol . l2012bh and smooth particle (SPH; 
Pakmor et al.ll2008l : Eiu et alJboiaY simulations, these non-degenerate companions should survive 
the SN impact and could be detectable. Thus, one of the simplest ways to distinguish between the 
SDS and DDS is to search for post-impact remnant stars (PIRSs) in Type la supernova remnants 
(la SNRs). 



Recent PIRS searches have studied two Galactic la SNRs fSN 1572. | Ruiz-Lapuente et a 



Ihara et al-lboOTt 



Kerzendorf et al 



200 



9: 



Gonza ez Hernandez et al 



Gonzalez Hernandez et al 



2012; Kerzendorf et al 



Edwards et al.l 2012 



200 



9; 



Kerzendorf et al 



2004; 



2012b . 



2012ah and tw o la SNRs in the 



Schaefer fc Pagnottal 12012 ). So far, only M 



and SN 1006, 

Large Magellanic Cloud (LMC; 
dwarf stars and the subgiant Tycho G star have emerged as possible PIRS candidates, and they are 
not well matched with companion models in the standard SDS channels, a result that may favor 
the DPS. However, the proper ties of a PIRS could ch ange significantly after the SN impact. For 
instance, Marietta et al .1 ()200Cl ) and IPan et al.l (|2012bl ) have shown that almost all the envelope of 
the RG in the RG-WD channel should be removed during the SN impa ct, and ~ 10 — 20% o f the 
MS star mass should be stripped and ablated in the MS- WD channel (Marietta et aliboool: Pan 
et al. l2012bl: \Uu et ali boid). Therefore, the PIRS in the RG-WD channel could be just a helium 
degenerate or non-degenerate core star with little hydrogen-rich envelope, where the former would 



-3 - 



be a pre- WD and the latter an extreme horizontal branch star. We note that the SN impact will 
not only strip and ab late the mass of a non-degenerate companion star, but also compress and 
deposit energy into it (|Pan et al.ll2012bl ). 



The evolu tion of PIRSs has bee n studied by iPodsiadlowskil (J2003J) for a 1M Q sub-giant com- 
panion and by IShappee et all (12013al ) for a 1M© MS companion. They found that PIRSs could be 
over-luminous due to the energy release from the SN energy deposition, suggesting that the SDS 
should be ruled out for several la SNRs. However, in their calculations, these authors assumed 
a fixed amount of energy input and mass stripping without performing detailed hydrodynamical 
calculations and, therefore, did not accurately ca lculate the shock compression in the stellar in- 
terior and the depth of the energy deposition. In |Pan et al.l (|2012al ). we studied the evolution of 
PIRSs using a detailed treatment of SN impact via three-dimensional hydrodynamics simulations. 
These simulations included the symmetry-breaking effects of orbital motion, rotation of the non- 
degenerate companions, and Roche-lobe overflow for the MS- WD channel. These three-dimensional 
simulation results were mapped into a one-dimensional stellar evolution code to simulate the post- 
impact evolution. It was found that MS-like PIRSs evolve to become subgiants (L ~ 10 — lOOL©) 
after a few hundred years and could be slowly rotating after stellar expansion. Although the model 
closest to the Tycho G star in these calculations was twice as bright, these results provide some 
support for Tycho G as a possible PIRS in the SDS. 

A new subcl ass of sub-luminou s SNe la, namely Type lax supernovae (SNe lax), recently has 
been proposed by lFoley et al.l (12012 ). This population could originate from the He- WD channel in 
the S DS via a helium double-detonation explosion or bv merger of a He WD with a CO WD ( Foley 



20121 : IWang et al.l 120131 )." The He- WD channel naturally explains the absence of hydrogen 



et al. 

lines, and two SNe lax have shown helium line s in their spe ctra, suggesting that helium must be 
present in the progenitor systems. iPan et al.l (|2010l . l2012bl ) have shown that only about < 5% 
of the mass of the helium star is lost into the SNR in the He- WD channel, an amount that is 
much lower than the hydrogen mass lost in the MS-WD channel and the RG-WD channel. The 
He- WD channel mainly contribute s to the prompt part of the delay-time distribution in population 
synthesis studies (~ 45 — 220 Myr, IWang fc Hanll2010bl ). and this is consistent with the distribution 
of SNe lax, since no SNe lax have been observed in elliptical galaxies. Finally, several He-WD 
binary systems wi th the properties required to be SN la p rogenitors have been observed, for ex ample 



200 



Geier et al 



200 



7n. CD-30 11223 dVennes et alJboi j_ Oeier 



KPD 1930+2752 (IMaxted et al 

et al. I2OI2I), and RX J0648.0-4418 ( |Mereghetti et ali bood). These systems have either a massive 
WD (M WD > 1.3M ) or a short orbital peri od (P mh < 0-05 d ). The helium nova V445 Pup is also 



considered to be a likely progenitor system (jKato et al.ll2008l ) 



In this paper, we follow the methods in IPan et al.l (|2012al ) and extend the study of post- 
impact evolution to helium-rich PIRSs (He PIRSs) in the SDS. We show that He PIRSs become 
hot blue-sub dwarf-like (sdO-like) stars after the SN impact and should be detectable in la SNRs. 
Hot blue subdwarfs are core helium-burning stars, and about half of t hem are in b inary systems. 
The formation of single hot blue sub-dwarfs is still an open question (jHeberl I2Q09T ) . and the He- 



- 4 - 



WD channe l could be one sourc e. He PIRSs may also contribute to the hypervelocity star (HVS) 



population (|Wang Sz Hanll2010d ). since the original orbital speed is very high at the time of the SN 



in our Galaxy; for example, US 708 has usini = 708 ± 15 km s 

consid er four different helium star progenitor models from the He- WD channel in I Wang et al 



(Hirsch et al. 


2005). We 


3 channel in 


Wane et al. 



420091) (HeWDa, HeWDb, HeWDc, and HeWDd in IPan et all (l2010f ). hereafter Paper I) and 
perform three-dime nsional hydrodyna mics simulations of the impact of SN la ejecta using the 
method described in lPan et al.l (|2012bl ) (hereafter Paper II). Based on the hydrodynamics results, 
w e carry out post-im pact stellar evolution simulations using a method similar to the one described 



m 



Pan et al.l (]2012al ) (hereafter Paper III). In the next section, the numerical codes and methods 
are described. Section [3] gives a detailed description of the He star progenitor systems in our 
calculations. In Section [H we present the simulation results and the possible evolutionary tracks of 
He PIRSs. We discuss several effects on the post-impact evolution and the possibility of He PIRSs 
as sources of single hot blue subdwarfs and/or HVSs in Section[5l In the final section, we summarize 
our simulation results and conclude. 



2. NUMERICAL CODES AND METHODS 



The simulation codes used in this paper are essentially the same a s those in Papers I I & III, 
inclu ding the three-dimensional hydrodynamics code FLASlfl version 3 (IFrvxell et al.ll2000l: Dubev 
et al. 



20081 ). to simulate the impact of SN la ejecta on the binary compan ions, and the stellar 



evolu tion code, MESaEI (Modules for Experiments in Stellar Astrophysics; iPaxton et al.l l201ll 



20131 ). to create the progenitor models and simulate the post-impact evolution. We perform the SN 
la explosion simulations and post-impact stellar evolution calculations using a technique similar 
to the one described in Paper II, but we focus on the He- WD channel. To link FLASH'S output 
with MESA's initial stellar models, we solve the hydrostatic equilibrium equation using a fourth- 
order Runge-Kutta solver with adaptive stepsize control, as described in Paper III. However, the 
assumption of ideal gas plus radiation equation of state (EOS) is no longer valid for He star 
companions, where the central density and temperature are much higher than in main sequence- 
like stars. We therefore update the EOS solver to include the OPAL, SCVH, and HELM EOS 
tables from the eos and kap modules in MESA. 



http : //f lash.uchicago . edu| 
^jhttp: //mesa, sourcef orge .net] 
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3. HE STAR PROGENITOR SYSTEMS 

The progenitor models are taken from the HeWDa, HeWDb, HeWDc, and HeWDd helium 
star models in Paper I, but without the simplification of uniform composition. These four models 
were generated with initial masses equal to 1.25, 1.35, 1.4, and 1.8 Mq and initial metallicity 
Z = 0.02. An artificial constant mass loss rate was adopted such that the evolutio n times and 



final h elium-star masses were consistent with the detailed binary evolution models of I Wang et al 



(|2009i ). Once the mass-loss phase ended, the stellar models were taken as initial models for the 
three-dimensional FLASH simulations. The physical properties of these four initial models are 
summarized in Table [U All four helium star models were relaxed on the three-dimensional grid by 
artificially damping the gas velocity for five dynamical timescales, ensuring that our models started 
in hydrostatic equilibrium and reducing the geometrical distortion introduced by passing from one 
dimension to three dimensions. The entropy, composition, density, and temperature of these four 
helium stars at the onset of the SN la explosion are shown in Figure [TJ 

Using the numerical setup and initial conditions for SN la explosions described in Paper II, we 
performed three-dimensional FLASH simulations of SN la explosions in close binary systems with 
resolutions of 6/8 AMR levels (equivalent to a 1024 3 uniform grid; see Paper II for definition). The 
initial binary syste ms were assumed to b e in Roche-lobe overflow (RLOF), and the SN model used 
was the W7 model (jNomoto et al.lll984l ). Although SNe lax are mostly sub-Chandrasekhar mass 
explosions, we chose the standard W7 explosion as the first case for comparison with our previous 
studies in Papers I and II. Since the mass loss from the non-degenerate star is sensitive to the 
detailed numerical setup of the SN model, we adopted the "W7 SN" sub-grid model, which has a 
power-law density distribution and a constant temperature distribution in radius that matches the 
exploding mass and energy in the W7 model (see Paper II for detailed descriptions). 

After the SN la impact, the helium star binary companion loses about ~ 4 — 6% of its mass 
and is heated, the degree to which depends on the progenitor model (see Table [2]). In addition, 
because of the shock interaction during the SN la impact, the central density of the helium star 

Table 1: The progenitor models at the onset of the SN explosion 



Model 


M (M ) 


t R (10 10 cm) 


lo gl0 L (Lq) 


logi T cff:0 (K) 


-Porb.o (sec) 


HeWDa 


0.697 


0.63 


1.30 


4.61 


1,033 


HeWDb 


0.803 


1.10 


1.62 


4.57 


2,257 


HeWDc 


1.007 


1.35 


2.11 


4.64 


2,682 


HeWDd 


1.206 


1.61 


1.09 


4.35 


3,410 



^The mass (Mo), radius (Ro), luminosity (Lo), effective temperature (T e ff,o), and or bital period (Porb.o) for different 
He star progenitor models at the time of the SN explosion, using the final masses in Wang et al. 1 20091 ). 
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Fig. 1. — The initial specific entropy, helium composition (Y), density, and temperature profiles as 
functions of the fractional mass before the SN la explosion for models HeWDa, HeWDb, HeWDc, 
and HeWDd in Table E 
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decreases to ~ 20 — 30% of its original value, and the central temperature decreases by ~ 10 — 15%, 
except for model HeWDa. Before the SN impact, the HeWDa model has a lower central helium 
composition than the other three models due to central helium burning, causing a higher positive 
entropy gradient at about m/M* ~ 0.5. Therefore, in this model the central density decreases 22% 
during the SN impact, while the central temperature increases by 12%. Figure [2] shows the angle- 
averaged one-dimensional radial post-impact profiles of the specific entropy, helium composition, 
density and temperature at the end of the FLASH simulations. 

With the eos (equation of state) and kap (opacity) module in MESA, we used the post-impact 
specific entropy profiles and composition profiles to construct hydrostatic models by solving the 
continuity and hydrostatic equations for the density p, pressure P, and radius r as functions of the 
enclosed mass m: 

— = — (1) 

dm 4irr 2 p 

dP Gm 



dm 4-7rr 4 

The entropy of each mass element was kept fixed, except that the entropy profiles were flattened 
in the outermost region (0.995 < m/M* < 1) to avoid negative entropy gradients. For the HeWDa 
model, the composition profile was adjusted to a uniform distribution due to the strong mixing 
during the SN impact and to avoid negative entropy and composition gradients. The hydrostatic 
solutions were taken as initial conditions for the models used in MESA. 

The initial luminosity profile for MESA was estimated using the radiative temperature gradient 
expression, 

where k is the opacity, a is the radiation constant, and c is the speed of light. Since our initial 
luminosity profile is based on an assumption of radiative equilibrium and our surface profiles are 
not as sharp as the standard surface profile in MESA, the calculated photospheric luminosity and 
effective temperature were very rough at the beginning and needed to be relaxed in MESA. A fixed 
time step, At = 10 -8 years, was enforced for the first 10 -6 years to relax the models. After 10 -6 
years, we allowed the time step to be automatically determined in MESA. 

Figure [3] shows the specific entropy, helium composition, density, and temperature profiles 
of the relaxed helium star models in MESA. The relaxed hydrostatic helium star models differ 
somewhat in the photospheric luminosity and effective temperature, but the stellar radius and 
interior density and temperature profiles are nearly the same as in the original models. Figure H] 
shows the changes of mass and radius before and after SN impact. Although He PIRSs only lose 
a few percent of their masses, post-impact radii increase by a factor of ~ 4 (see Table [2|). These 
changes in radius dramatically alter the He PIRSs. 
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Fig. 2. — Similar to Figure Q] but for post-impact angle-averaged profiles at the end of the FLASH 
simulations. 
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Fig. 3. — Relaxed post-impact companion models in MESA for specific entropy, helium composition 
(V), density, temperature, and enclosed luminosity as functions of the fractional mass for models 
HeWDa, HeWDb, HeWDc, and HeWDd in Table 
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Fig. 4. — Stellar mass and radius for all considered helium star models before (crosses) and after 
(filled squares) supernova impact. 



4. POST-IMPACT EVOLUTION 



Once we have the relaxed hydrostatic stellar models, the post-impact evolution of He PIRSs 
can be easily calculated in MESA. In this section, we describe this evolution. Hertzsprung-Russell 
(H-R) diagrams of post-impact evolutionary tracks of He PIRSs are presented in luminosity versus 
effective temperature, surface gravity versus effective temperature, and color-magnitude forms, 
providing diagnostics for searches in future observations for He PIRSs in la SNRs. 



4.1. Evolutionary Tracks 

Similarly to the MS-like PIRSs in Paper III, the He PIRSs expand rapidly and dramatically 
due to the release of energy deposited by the SN impact. Figure [5] shows the evolution of the 
photospheric radius, luminosity, and effective temperature as functions of time. All He PIRSs 
expand on a timescale of ~ 10 — 30 years, depending on their pr ogenitor models. These ex pansion 



rates are determined by the local radiative diffusion timescale (jHenvev &; LT^cuverlll969l ). which 
is associated with not only the stellar structure, but also the amount and depth of SN energy 
deposition. Since He stars are more compact than MS-like stars, the depths of energy deposition 
are shallower, causing a shorter local radiative diffusion timescale. Therefore, heat transfer initially 
occurs more rapidly than the thermal expansion. The effective temperature starts to increase at 
~ 10~ 2 — lO^ 1 years and then continues increasing for another 10° — 10 1 yr up to ~ 30, 000—50, 000 K. 
Subsequently, the surfaces of the He PIRSs cool off to ~ 10, 000 — 30, 000 K on a timescale of 
~ 10 years due to the expansion. The stars quickly become luminous helium OB stars {L ~ 
10 3 -10 4 L Q ). When the deposited energy has radiated away, the photospheric radii reach maximum 
values at around ~ 10 — 30 yrs. Subsequently, the stars contract and release gravitational energy, 
turning them into sdO-like stars for ~ 10 7 years. After ~ 10 7 years, the core helium of He PIRSs 
should be exhausted, and the stars should evolve to the helium red giant phase. The simulations 
are terminated at ~ 10 5 years, since the SN la remnant may not be recognizable after this time. 

We have noticed that there is < 1% artificial mass loss during angle-averaging of the post- 
impact density, since the helium star is close to the edge of the simulation box at the end of the 
FLASH simulations. However, the most important factors that control the post-impact evolution 
are the amount of energy deposited and the corresponding depth, which are mainly concentrated 
in the outer 5% of mass (see the temperature bump in Figure [2]) . We have tested the effect of mass 
loss by changing the mass of the post-impact remnant star and conclude that the artificial mass 
loss will not lead to a notable difference in the post-impact evolution. 

The evolutionary tracks of He PIRSs in the H-R diagram are plotted in Figures [6] and [7J 
representing the effective temperature versus luminosity and effective temperature versus surface 
gravity, respectively. Based on our post-impact simulations, we predict that He PIRSs will evolve 
to luminous OB stars on a timescale of ~ 30 years and will fade within a hundred years. Therefore, 
only in young SN la remnants can one observe luminous helium OB stars at the center of the SNR 
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Fig. 5. — Evolution of the photospheric radius, effective temperature, and luminosity as functions 
of time. Each line shows the evolution of a He PIRS in Table [2j 
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if the non-degenerate companion was a He star. However, an sdO-like star is observable for the 
remai ning time. Note that for the He- WD channel in the SDS the delay time (~ 45—220 Mvr. Wang 



& Han l2010bl ) is much shorter than in other SDS channels since the helium star was formed from a 
more massive star. This suggests that He PIRSs are expected to be detected only in star-forming 
regions. 

As SNe lax represent a sub-class of sub-luminous SNe la, their explosion energy and ejecta 
speed are lower than in the standard W7 model, causing a reduced effect of the SN impact on 
binary companions. In Paper III, we studied the effect of the SN explosion energy on the post- 
impact evolution and found that the explosion energy dramatically affects the amount and depth of 
SN energy depositition in MS-like PIRSs. For a lower explosion energy, the effect of the SN impact 
is shallower, causing a shorter radiative diffusion timescale. However, these differences become 
small once the deposited energy has radiated away. Therefore, for He PIRSs with sub-luminous 
explosions, He PIRSs behave similarly after ~ 30 years, but with less mass lost during the SN 
impact. 



4.2. Color- Magnitude Diagram 

For direct comparison with observations, we convert the luminosity to magnitude in the optical 
bands. Given the effective temperature (T c g) and photospheric radius (R) of a PIRS, the magnitude 
of the PIRS can be calculated using Equation [J] with an assumption of blackbody radiation: 



m S, 



-2.5 log 



10 



J S x (irB x )d\ (R 
J(fSc/\*)S x d\ \d 



(4) 



where S\ is the sensitivity function of a given filter at wavelength A, B\ is the Planck function, d is 
the distance of the PIRS, and /° = 3.631 x 10 -20 erg cm" 2 s" 1 Hz" 1 is the zero-point value in the 
AB magnitude system. Figure [5] shows the color-magnitude evolutionary trajectories of He PIRSs. 
The absolute magnitudes are calculated using the broadband u and v filters. 



5. DISCUSSION 

In this section, we study the possible observational effects of the post SN impact evolution. 
In particular, we examine the effect of nickel contaminations from the supernova ejecta on the 
post-impact evolution. The surface rotational speeds of He PIRSs during post-impact evolution 
are also predicted as functions of time. Furthermore, we also discuss the possibility that He PIRSs 
could be a source of single sdO stars and/or HVSs. 
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Fig. 6. — Evolutionary tracks in the H-R diagram for different He PIRS models. Each line represents 
an evolutionary track of a He PIRS in Table [2] over an interval of 10 5 years. The filled squares 
indicate the conditions at 0.1 years after the SN impact; filled stars, 10 years after the SN impact; 
filled triangles, 10 2 years after the SN impact; filled circles, 10 3 years after the SN impact. The 
He PIRS models with letter "n" represent the cases without nickel contamination (with dashed 
lines) . 
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Table 2: The progenitor models immediately after the SN impact 



Model 


M SN t 


(AM /Mo) 




(Rsn/Ro) 


l°g 10 -^SN 


loglO ^eff.SN 


^linear.SN 


M Ni 




(M Q ) 




(10 10 cm) 




(L & ) 


(K) 


(km sec -1 ) 


(1O~ 4 M ) 


HeWDa 


0.656 


5.88 % 


2.71 


4.30 


0.265 


4.03 


734 


15.0 


HeWDb 


0.748 


6.85 % 


3.93 


3.58 


0.510 


4.01 


550 


2.38 


HeWDc 


0.962 


4.47 % 


5.63 


4.17 


0.792 


4.01 


509 


5.64 


HeWDd 


1.126 


6.6 3% 


7.24 


4.50 


1.16 


4.03 


446 


1.75 



The mass (Msn), mass change (AM = Msn — Mo ), radius (Rsn), luminosity (Lsn), effective temperature (Teff.SN), 
linear spatial velocity (i>ii I1C ar.SN), and mass of bound nickel (MNi) of initial relaxed post-impact hydrostatic models 
in MESA. 
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Fig. 7. — Similar to Figure [6] but for surface gravity (in cgs units) versus effective temperature (K). 
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Fig. 8. — Similar to Figure [6] but for V magnitude versus B-V color. The magnitudes are absolute 
magnitudes using the broadband u and v filters in the AB magnitude system. 
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5.1. Nickel Contamination 

The envelope of a companion star could be contaminated by the supernova ejecta during the 
supernova impact or as fallback. In Paper II, we have shown that this nickel contamination is 
much greater for He PIRSs (1.7 x 1O" 4 M - 1.5 x 1O~ 3 M ; see Table ED than MS-like PIRSs 
(< 10~ 5 Mq). We note that the nickel contamination can affect the post-impact evolution for 
He PIRSs, since the change of metallicity in the envelope will also change the opacity and radiative 
diffusion timescale. To study this effect, two He PIRSs (labeled "HeWDbn" and "HeWDcn") 
without nickel contamination have been reconstructed by removing the bound nickel at the end 
of the FLASH simulations. The hydrostatic profiles of the reconstructed models do not change 
significantly, but they have slightly smaller radii. Figure [9] shows the evolution of the photospheric 
radius, effective temperature, and luminosity of these two models with and without (cases with 
the letter "n") nickel contamination. Removing the bound nickel causes the opacity to be lower 
in the outer regions, making stars more transparent and causing smaller photospheric radii. It is 
clear that the maximum radii of He PIRSs without nickel contamination are much smaller than 
the same He PIRSs with nickel contamination. Note that we used the fixed metal tables (Type I) 
opacity module in MESA. Thus, the opacity contributed by heavy metals in our calculations is not 
very accurate, since the abundance ratios are different in the nickel-contaminated region than the 
abundances assumed in OPAL. However, the difference between models with and without nickel 
contamination is not very significant in the H-R diagrams (see Figured]). Thus, these differences 
can be treated as an uncertainty in our post-impact simulations. 

Energy generation by 28 Ni decay and 28 Co decay in the SN ejecta-contaminated regions is 
neglected in our calculations. This is justified by the fact that the energy generated by these 
nuclear decays is much smaller than the energy deposited from the supernova ejecta since the mass 
of bound nickel is small. 



5.2. Surface Rotational Speed 

In Paper III, we showed that the surface rotational speed could decrease to ~ 10 — 20 km s _1 
for MS-like PIRSs, if the specific angular momentum is conserved during the post-impact evolution. 
This result implied that the possible PIRS candidate Tycho G could not be completely ruled out 
because of its low surface rotation speed. Here, we apply the same method to He PIRSs. The 
He stars were set into rigid rotation at the beginning of the FLASH simulations becaus e the 



synch ronization time due to tidal locking for a given mass ratio q and orbital period P (IZahn 
19771 ). tgync ~ 10 4 ((1 + q)/2q) 2 P 4 yrs <C 1 day, is extremely short in the He- WD channel. After 



the SN impact, a He PIRS loses < 10% of its angular momentum and is no longer in a state of 
rigid-body rotation. Thus, an angle-averaged, post-impact radial angular velocity profile can be 
calculated by averaging different latitudes and longitudes in spherical coordinates in the FLASH 
output. The post-impact surface rotational speed can be calculated by assuming the specific angular 
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Fig. 9. — Similar to Figure but for HeWDb and HeWDc with and without nickel contaminations. 
The letter "n" in He PIRS models represents the cases without nickel contamination. 
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momentum to be conserved after the SN impact. Figure [TOl shows the post-impact evolution of the 
surface rotational speed versus evolution time. It is found that the angle-averaged angular velocity 
profiles exhibit some variations in the envelope region, and we estimate the uncertainty by using 
the standard deviation of the specific angular momentum within the envelope. Since He stars are 
more compact and are closer to the accreting WD at the time of the SN explosion, the surface 
rotational speed should be much higher than for MS-like companions in the MS- WD channel. The 
surface rotational speed at the time of the SN la explosion could be as high as ~ 200 km s- 1 if 
the He star companion is close to co-rotation due to tidal locking. For He PIRSs, our calculations 
show that the surface rotational speed should decrease to < 10km s~ 1 when the star expands to its 
maximum size. However, because the radiative diffusion timescale is short, the surface rotational 
speed is slow only during the first hundred years and will eventually increase to ~ 100 km s~ 1 after 
a few hundred years. Therefore, if He PIRSs exist in historical la SNRs, their surface rotational 
speeds could be higher than 50 km s _1 , depending on the progenitor models and the ages of the 
SNR. 



5.3. He PIRSs as hypervelocity stars? 



During the past decade a number of HVSs have been observed in the halo of our Galaxy. 
They are sub-luminous O- or B-type stars or MS B stars with high radial velocities that could 
exceed the escape velocity of the Milky Way. The common explanation for their formation involves 
tidal ejectio n in binary stars (or triple stars) a ssociated with ei ther a massive black hole or binary 
black holes (jHillsl Il988l : lYu fc Tremainel hood ) . Alternatively, Iwang fc Han] (|2009l ) suggest that 
the surviving He PIRSs in SNe la within the He- WD channel could be a source of HVSs. In our 
calculations, He PIRSs are sdO-like stars for ~ 10 7 yr and thereafter evolve to the helium red-giant 
phase. The final linear velocity (original orbital velocity plus kick velocity) of our He PIRSs is 
~ 400 — 800 km s _1 (see Table [2]). Assuming the la SNRs will not be recognizable after 10 5 years, 
He PIRSs could move to a distance of ~ 10 kpc within 10 7 years, becoming single sdO-like stars in 
the halo of our Galaxy. The HVS US 708 (or HVS2) has a radial velocity v sini = 708 ± 15 km s _1 , 
T c fj = 44. 500 K. log a = 5.23. and Galactic latitude b = +47.05° at a distance of 19 kpc ( Hirsch 
et al. 



200 



5j), giving a displacement of 14 kpc from the Galactic plane. The effective temperature 



and surface gravity of US 708 are consistent with our He PIRS models, but they require a subdwarf 
lifetime longer than several times 10 7 yr to reach such a distance, if the SN exploded in the Galactic 
plane. Therefore, a less massive model may better match US 708, since less massive models have 
higher linear speeds and longer subdwarf lifetimes. 



6. CONCLUSIONS 



We have investigated the post-impact evolutio n of He PIRSs wit hin the SDS for SNe la via 
numerical simulations. Four helium star models from I Wang et al. ( 2009 ) are considered in our calcu- 
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Fig. 10. — Evolution of surface rotational speed for all He- WD models in Table El The error is 
estimated using the variation of specific angular momentum within the envelope of initial hydro- 
static solutions (see Section [5.21 for detailed description). The number of error bars is only chosen 
for visualization and does not represent the actual number of data points. 
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lations. We performed three-dimensional hydrodynamics simulations using the methods that were 
described in Paper II. The post-impact evolution of these stars has been studied by reconstruct- 
ing hydrostatic models based upon the final output in the hydrodynamics simulations and then 
interpolating them into a one-dimensional stellar evolution code, MESA. It is found that He PIRSs 
expand dramatically and evolve to become luminous OB stars within about ~ 10 — 30 years after 
the SN la explosion. This phase is short (< 100 years), and therefore these luminous OB stars are 
not likely to be detected in historical la SNRs. After ~ 30 years, He PIRSs contract and evolve into 
hot blue-subdwarf-like (sdO-like) stars due to the release of gravitational energy. Therefore, we 
predict that most He PIRSs should be sdO-like stars and could be detectable in nearby la SNRs. 

We also predict that these He PIRSs should be rapidly rotating (v rri t 50 km s" 1 ). Al though 



a few fast rotating hot blue-subdwarfs have been observed rece ntly (IGeier et a 



single hot blue-subdwarfs (sdOs and sdBs) are slowly rotating (IGeier &: Heber 



Hebei 



2011 



20131 ). most 



2013). If single hot 
20091 ) . merger via a 



blue-subdwarfs originate from the merger of two He WDs (sdOs/sdBs ; 
common envelope phase on the red giant branch (jPolitano et al.ll2008l ). or the He- WD channel in 
SNe la (sdOs), the theoretical predictions cannot explain these observations, suggesting that other 
mechanisms operate to slow down the rotation, such as magnetic braking. 

The He- WD binary channel is favored for the prompt DTD in the SDS and is expected to occur 
in star-forming regions. Since the orbital period immediately prior to the SN la explosion in the 
He- WD channel is extremely short (< 1 hour), the system is expected to be tidally locked. Thus, 
He stars should be rapidly rotating at the time of the SN la explosion. Although some angular 
momentum will be lost during the SN impact, He PIRSs are still expected to be rapidly rotating 
after ~ 30 years. The spatial velocity of He PIRSs is also expected to be high, reflecting the high 
orbital speed at the time of the SN la explosion and implying that He PIRSs could contribute to 
the HVS population (i.e. US 708). 



The simulations presented here were carried out using the NSF XSEDE Ranger system at 
the Texas Advanced Computing Center under allocation TG-AST040034N. FLASH was developed 
largely by the DOE-supported ASC/Alliances Center for Astrophysical Thermonuclear Flashes 
at the University of Chicago. This work was partially supported by the Computational Science 
and Engineering (CSE) fellowship at the University of Illinois at Urb ana-Champaign. Part of the 



analysis of FLASH data was completed using the analysis toolkit yt ([Turk et al 



2011 



)• 
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